UCT-TP-279/2010 
February 2010 



o 



The Operator Product Expansion Beyond Perturbation Theory 

in QCD 1 







O 

(N 

On ' C. A. Dominguez 

^ ! 

' Centre for Theoretical Physics and Astrophysics 

- ■ University of Cape Town, Rondebosch 7700, South Africa, and Department of Physics, 



SteUenbosch University, Stellenbosch 7600, South Africa 



Abstract 

> 

I The Operator Product Expansion (OPE) of current correlators at short distances beyond per- 

' turbation theory in QCD, together with Cauchy's theorem in the complex energy plane, are the 

fSJ ■ pillars of the method of QCD sum rules. This technique provides an analytic tool to relate QCD 

I with hadronic physics at low and intermediate energies. It has been in use for over thirty years 

to determine hadronic parameters, form factors, and QCD parameters such as the quark masses, 
and the running strong coupling at the scale of the r-lepton. QCD sum rules provide a powerful 
complement to numerical simulations of QCD on the lattice. In this talk a short review of the 
^ ' method is presented for non experts, followed by three examples of recent applications. 
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1 Introduction 



The method of QCD sum rules, introduced by Shifman, Vainshtein, and Zakharov [Ij more than 
thirty years ago, has become a powerful technique to study hadronic physics in the low energy 
resonance region by means of QCD [2] . It is also a complementary tool to numerical simulations 
of QCD on a lattice. The range of applications has steadily grown over the years and covers 
the determination of the hadronic spectrum (masses, couplings, and widths), electromagnetic, 
weak, and strong form factors, quark masses and the strong coupling, the extension of QCD 
to finite temperature and its related topics of chiral symmetry restoration and quark-gluon 
deconfinement [3]. It is also a tool to confront QCD predictions with experimental data, e.g. 
from e~^e~ annihilation and hadronic r-lepton decays. This method is based on two fundamental 
pillars: (i) the operator product expansion of current correlators at short distances, extended 
beyond perturbation theory to incorporate quark-gluon confinement, and (ii) Cauchy's theorem 
in the complex energy (squared) plane, often referred to as quark-hadron duality. To be more 
specific, let us consider a typical object in QCD in the form of the two-point function, or current 
correlator 



where the local current J{x) is built from the quark and gluon fields entering the QCD La- 
grangian, and it has definite quantum numbers. Equivalently, this current can be written in 
terms of hadronic fields having the same quantum numbers. The specific choice of current will 
depend on the application one has in mind. For instance, if one is interested in determining the 
hadronic properties of the /?"'"-meson, then one would choose the QCD vector isovector current 
J^{x) = d(x)7^n(x), and its hadronic realization in terms of the /9+-meson field. If the goal 
is to determine the values of the light quark masses, then the ideal object would be the corre- 
lator involving the axial-vector current divergences J5(x)|*- = {nii + mj)-ip^{x)j5ipj{x), with i,j 
the up, down, or strange quark flavors. The hadronic representation of this correlator contains 
the pseudoscalar meson (vr or K) mass and coupling, followed by its radial excitations and the 
hadronic continuum. The tool to relate these two representations is Cauchy's theorem in the 
complex energy (squared) plane, to be discussed shortly. 

The QCD correlator, Eq.(l), will contain a perturbative piece (PQCD), computed up to a given 
loop order in perturbation theory, and a non perturbative part mostly reflecting quark-gluon 
confinement. The leading order in PQCD is shown in Fig.l. Since QCD has never been solved 
analytically, the effects due to confinement can only be introduced by parameterizing quark and 
gluon propagator corrections effectively in terms of vacuum condensates. This is done as follows. 
In the case of the quark propagator 



the quark propagator correction S(p^) would contain the information on confinement. One 
expects this correction to peak at and near the quark mass-shell, i.e. for p ~ in the case of 
light quarks. Effectively, this can be viewed as in Fig. 2, where the (infrared) quarks in the 
loop have zero momentum and interact strongly with the physical QCD vacuum. This effect is 
then parameterized in terms of the quark condensate {0\q{0)q{0)\0) . Similarly, in the case of the 
gluon propagator one would have 
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Figure 1: Leading order PQCD correlator. All values of the four-momentum of the quark in the 
loop are allowed. The blue wiggly line represents the current of momentum q {—q^ >> 0). 



where the gluon propagator correction will peak at A; ~ 0, and the effect of confinement in this 
case will be parameterized by the gluon condensate (0|as G^*^ • G^y|0) (see Fig.2). In addition 



Figure 2: Quark propagator modification due to (infrared) quarks interacting with the physical 
QCD vacuum, and involving the quark condensate. Large momentum flows through the bottom 
propagator. 



to the quark and the gluon condensate there is a plethora of higher order condensates entering 
the OPE of the current correlator at short distances, i.e. 

W)\qCD = Coi C2N+2{q\l^^) (0|d2iV+2(;u')|0) , (4) 

where /i^ is the renormalization scale, and where the Wilson coefficients in this expansion depend 
on the Lorentz indices and quantum numbers of J{x) and of the local gauge invariant operators 
On built from the quark and gluon fields. These operators are ordered by increasing dimension- 
ality and the Wilson coefficients, calculable in PQCD, fall off by corresponding powers of —(f'- 
Since there are no gauge invariant operators of dimension d = 2 involving the quark and gluon 
fields in QCD, it is normally assumed that the OPE starts at dimension d = 4 (with the quark 
condensate being multiplied by the quark mass). This is supported by results from QCD sum 
rule analyses of r-lepton decay data, which show no evidence of d = 2 operators [1]. The unit 




Figure 3: Gluon propagator modification due to (infrared) gluons interacting with the physical 
QCD vacuum, and involving the gluon condensate. Large momentum flows through the quark 
propagators. 

operator in Eq.(4) has dimension d = and CqI stands for the purely perturbative contribution. 
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The Wilson coefficients as well as the vacuum condensates depend on the renormalization scale. 
In the case of the leading d = 4 terms in Eq.(4) the /i^ dependence of the quark mass cancels 
the corresponding dependence of the quark condensate, so that this contribution is a renormal- 
ization group (RG) invariant. Similarly, the gluon condensate is also a RG invariant, hence once 
determined in some channel these condensates can be used throughout. At dimension d = 6 
there appears the four-quark condensate, obtained from Fig.l at the next to leading order (one 
gluon exchange) and allowing all four quark lines to interact with the physical vacuum (see 
Fig. 4). While this condensate has a residual renormalization scale dependence, this is so small 
that in practice it can be ignored. The four-quark condensate, while relatively small, is crucial 
to explain the large p{770) - ai(1260) mass splitting. In most applications —q^ is chosen large 
enough so that the condensates of higher dimension (d > 8) can be safely ignored. The numerical 




Figure 4: The four-quark condensate of dimension d = 6 in the OPE. This is responsible for the 
p — ai mass splitting. Large momentum flows through the gluon propagator. 



values of the vacuum condensates cannot be calculated analytically from first principles as this 
would be tantamount to solving QCD exactly. One exception is that of the quark condensate 
which enters in the Gell-Mann-Oakes-Renner relation, a QCD low energy theorem following 
from the global chiral symmetry of the QCD Lagrangian. Otherwise, it is possible to extract 
values for the leading vacuum condensates using QCD sum rules together with experimental 
data, e.g. e~^e~ annihilation into hadrons, and hadronic decays of the r-lepton. Alternatively, 
as lattice QCD improves in accuracy it should become a valuable source of information on these 
condensates. 

Turning to the hadronic sector, bound states and resonances appear in the complex energy 
(squared) plane (s-plane) as poles on the real axis, and singularities in the second Riemann 
sheet. In addition there will be multiple cuts reflecting non-resonant multi-particle production. 
All these singularities lead to a discontinuity across the positive real axis. Choosing an integra- 
tion contour as shown in Fig. 5, and given that there are no further singularities in the complex 
s-plane, Cauchy's theorem leads to the finite energy sum rule (FESR) 

[•so I If 

/ ds- f{s) ImU{s)\HAD = -— f ds f{s)U{s)\QCD , (5) 
Jsth 27^^ Jc(\so\) 

where /(s) is an arbitrary (analytic) function, sth is the hadronic threshold, and the finite radius 
of the circle, so, is large enough for QCD and the OPE to be used on the circle. Physical observ- 
ables determined from FESR should not depend on sq- In practice, though, this independence 
is not exact, and there is usually a region of stability where observables are fairly independent 
of typically in the range sq — 1 — 4 GeV^. The variation of an observable in the stability 
region is incorporated into the error of the determination. Equation (5) is the mathematical 
statement of what is usually referred to as quark-hadron duality. Since QCD is not valid in the 
time-like region (s > 0), in principle there is a possibility of problems on the circle near the real 
axis (duality violations). I shall come back to this issue later. The right hand side of this FESR 
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Figure 5: Integration contour in the complex s-plane. The discontinuity across the real axis 
brings in the hadronic spectral function, while integration around the circle involves the QCD 
correlator. 



involves the QCD correlator which is expressed in terms of the OPE as in Eq.(4). The left hand 
side calls for the hadronic spectral function which is written as 

ImIi{s)\HAD = ImU{s)\poLE + ImU{s)\REsd{so- s) + ImU{s)\pQCD 0{s - sq) , (6) 

where the ground state pole, absent in some channels, is followed by the resonances which merge 
smoothly into the hadronic continuum above some threshold sq. This continuum is expected 
to be well represented by PQCD if sq is large enough. Due to this, if one were to consider 
an integration contour in Eq.(5) extending to infinity, the cancellation between the hadronic 
continuum in the left hand side and the PQCD contribution to the right hand side, would 
render the sum rule a FESR. The performance of the contour integral in the complex s-plane is 
discussed in the next section. 



2 Finite energy sum rules 

The integration in the complex s-plane of the QCD correlator is usually carried out in two 
different ways. Fixed Order Perturbation Theory (FOPT) and Contour Improved Perturbation 
Theory (CIPT) [5]. The first method treats running quark masses and the strong coupling 
as fixed at a given value of sq. After integrating all logarithmic terms {ln(—s/fj?)) the RG 
improvement is achieved by setting the renormalization scale to /i'^ = —sq. In CIPT the RG im- 
provement is performed before integration, thus eliminating logarithmic terms, and the running 
quark masses and strong coupling are integrated (numerically) around the circle. This requires 
solving numerically the RGE for the quark masses and the coupling at each point on the circle. 
The FESR Eq.(5), with /(s) = 1, in FOPT can be written as 




M2N+2{so) 



(7) 
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where the dimensionless PQCD moments M2Ar_|_2(so) are given by 



M2N+2{so) 



/•So ]^ 

/ ds s'^ - ImU{s)\pQCD 
Jo 



(8) 



If the hadronic spectral function is known in some channel from experiment, then ImIl{s)\HAD = 
/mil (s) I DATA 5 and Eq.(7) can be used to determine the values of the vacuum condensates. 
Subsequently, Eq.(7) can be used in a different channel to determine the masses and couplings 
of the hadrons in that channel. It is important to mention that the correlator n(g^) is generally 
not a physical observable. However, this has no effect in FOPT as the unphysical constants in 
the correlator do not contribute to the integrals. The situation is quite different in GIFT where 
Eq.(5) cannot be used for unphysical correlators. For instance for a correlator whose physical 
counterpart is the second derivative (needed to eliminate a first degree polynomial), Cauchy's 
theorem and the resulting FESR must be written for the second derivative. In this case one has 
to use the following identity [6j 



where 



/ ds g{s) n(s) = / ds [F{s) - F{so)] n"(s) 

Jc{\so\) Jc{\so\) 



/s i-s' /-so 

ds' / ds"g{s") - / ds"g{s") 



(9) 



(10) 



and g{s) is an arbitrary function. This is easily proved by integrating by parts the right hand 
side of Eq.(9) and using Eq.(lO) to obtain the left hand side. In this case Eq.(5) becomes 

" ds g{s) - ImU{s)\HAD = - ^ i ds [F{s) - F{so)]I\" {s)\qcd , (11) 

2m Jc(\sa\) 



sth 



vr 



which is the master FESR to use in GIFT. The running quark masses and the running strong 
coupling entering 11" (s) are now functions of the integration variable and are not fixed as pre- 
viously in FOFT. The running coupling obeys the RGE 



s ^"''f = /3{as) = - /3n as{-s] 

N=0 



N+2 



ds 



(12) 



where = Os/tt, and e.g. for three quark flavors /3q = 9/4, /3i = 4, /32 = 3863/384, /S^, = 
(421797/54 + 3560C(3))/256. In the complex s-plane s = sq e"' with the angle x defined in the 
interval x € (— 7r,7r). The RGE then becomes 



das{x) 
dx 



Pn as (a 



\N+2 



(13) 



iV=0 



This RGE can be solved numerically at each point on the integration contour of Eq.(ll) using 
e.g. a modified Euler method, providing as input as{x = 0) = as{—so). Next, the RGE for the 
quark mass is given by 



s dm{—s) 
m ds 



E 

M=0 



IM % 



(14) 
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where e.g. for three quark flavors 70 = 1, 71 = 182/48, 72 = [8885/9 - 160C(3)]/64, 73 = 
[2977517/162 - 148720 C(3)/27 + 2160 C(4) - 8000C(5)/3]/256. With the aid of Eqs. (12)-(13) 
the above equation can be converted into a differential equation for m(x) and integrated, with 
the result ^ 

m{x) = m(0) exp| - if dx' 7m [as{x')]^'^+'^^ , (15) 

where the integration constant m(0) can be identified with m(so). 



3 APPLICATIONS 



As a first appHcation I sketch the determination of the strong couphng at the scale of the r 
mass [7J using experimental data on the ratio Rr = T{t — )• hadrons) /T{t — ?> leptons). The 
r-decays proceed through the V — A current, so that J{x) in Eq.(l) becomes J{x) — )■ J^{x) = 
d{x)j^{l — 75)n(x) in the non-strange sector. The ratio Rr can then be written as 



Rr = constant 



Mi 



So Sq 



(1 + 2—) Imn(^) (s) + Imn(°) (s) 
So 



(16) 



where H^'^) stands for the spin J = or J = 1 components. The QCD spectral functions are 
known in PQCD up to five-loop order, and the non-perturbative contribution is very tiny as it 
only starts at dimension d = 6 (due to the presence of the integration kernel). Hence, the strong 
coupling is basically the only unknown which can be determined after using the experimental 
data for the left hand side of the FESR above. The latest results are as{M'^) = 0.332 ± 0.016 
from [8] using FOPT and GIFT, and asiM"^) = 0.344 ± 0.009 from [9] using GIFT. Gomputing 
the evolution of the coupling to higher energies and comparing with independent determinations 
at the scale of the Z mass provides the most accurate test of asymptotic freedom. It also provides 
the most accurate and transparent determination of Ug at the five-loop level. 
Next I discuss an application where the integration kernel /(s) in Eq.(5) is of great importance 
[lOj . In the axial- vector channel, the FESR Eq.(7) with f{s) = 1 and = can be confronted 
with data from r-decay. The hadronic spectral function is then written as the sum of the pion 
pole and the resonance data known up to the kinematical end point sq = Af^. The moment 
M2(so) is known up to five- loop order in FQGD, so that the FESR can be used to confront 
the resonance data plus FQGD with e.g. f-,^. As seen from Fig. 6 the agreement is rather 
poor, except possibly near the end point. At first sight, this may be interpreted as a signal for 
quark-hadron duality violations near the real axis even at high enough energy. In fact, it has 
been known for quite some time that e.g. the Weinberg (chiral) sum rules are not saturated by 
the r decay data unless one introduces pinched integration kernels, e.g. /(s) = [1 — {s/sq)]^^^^^ 
|llj . Unfortunately, the r-lepton is not massive enough to probe higher energy regions. In 
spite of this it is still possible to enlarge the energy range by introducing as integration kernel 
a polynomial /(s) = P{s) designed to eliminate the (unknown) hadronic contribution to the 
integral between si and sq > si, where si is at or near the end point of the data. It has been 
shown |10j that the optimal degree of P{s) is the simplest, i.e. the linear function 



P{s) = 1 



2s 
So + si 



(17) 
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Figure 6: Results for from the standard FESR in the axial-vector channel, Eq. (7) with N=0, 
with no dimension d = 2 term, and using CIPT, with A = 365 MeV (as(M^) = 0.335). The 
straight line is the experimental value of /^r, and the points are the integrated data with the 
experimental errors. 



so that 



so 



constant x / P{s)ds = , 



(18) 



In this case the complete FESR becomes a linear combination of a dimension-two and a dimension- 
four FESR, which from Eqs.(7) and (17) it is given by 



2/^ 



ds P{s) - Imn{s)\DATA + -p^ 
Q TT 47r^ 



M2(so) 



2so 



So + Si 



M4(so) 



+ 



1 

4^ 



02(62) + 



So + Si 



^4(04) 



+ A(so), 



(19) 



where the pion pole has been separated from the data, and the chiral limit is understood. The 
term A(so) is the error being made by assuming that the data is constant in the interval si — sq. 
It is possible to estimate this error which turns out to be two to three orders of magnitude 
smaller than 2/^ on the left hand side of Eq.(19) [10]. As can be seen from Fig. 7 the FESR 
Eq.(19) shows an excellent consistency between QCD and the r data in the axial-vector channel 
in a remarkably wide region sq — 4 — 10 GeV^. A similar consistency is also found in the vector 
channel, where QCD is now confronted with zero (there is no pole in this channel). This result 
shows no evidence for quark-hadron duality violations in the vector and axial- vector channels. 
In addition, this method can prove quite useful in QCD FESR applications needing an extension 
of the energy region of available data. 

Finally, I discuss the use of FESR to determine the values of the QCD light quark masses [6]. 
As mentioned in the Introduction, in the case of the up-, down, and strange-quark masses the 
ideal current operator in Eq.(l) is the axial- vector current divergence. Its advantage is that the 
masses appear here as overall multiplicative factors, rather than as subleading power corrections 
like in other correlators, e.g. the vector or axial-vector correlators. The great disadvantage 
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Figure 7: Results for f.,^ from the FESR in the axial-vector channel, Eq. (19), with 02(02) = 0, 
C4(04) = 0.05 GeV^, , and using CIPT with A = 365 MeV {as{M^) = 0.335). The straight 
line is the experimental value of f-j^, and the points with errors are the integrated data up to 
si ~ M^. 



is that there is no direct experimental data beyond the pseudoscalar meson poles, i.e. the 
hadronic resonance spectral function, ImIl{s)\RES in Eq.(6) is not known experimentally. The 
only available information is that there are two radial excitations in the non-strange (vr) as 
well as in the strange {K) channel with known masses and widths. This is hardly enough 
to reconstruct the full spectral function. In fact, inelasticity, non-resonant background, and 
resonance interference are impossible to guess so that a model is needed for the resonant spectral 
function. This fact, which introduces a serious systematic uncertainty, has affected all quark 
mass determinations using QCD sum rules until recently [6]. The breakthrough has been to 
introduce an integration kernel in the FESR tuned to suppress substantially the resonance energy 
region above the ground state. This kernel is of the form /(s) = A5(s) = 1 — s — ai s^, where 
the coefficients are fixed by requiring that A5(s) vanish at the peak of the two radial excitations, 
i.e. A^{Mf) = A5(M|) = 0. This has the effect of reducing the resonance contribution to the 
FESR to a couple of a percent of the ground state contribution, well below the uncertainty due 
to Os- Hence the systematic uncertainty is essentially eliminated and one finds, using FOPT 
and CIPT 

m„(2 GeV) = 2.9 ± 0.2 MeV , 
md(2 GeV) = 5.3 ± 0.4 MeV , 

mud = ^ = 4.1 ± 0.2 MeV , 

ms{2 GeV) = 102 ± 8 MeV . 

This is at present the most accurate determination of the light quark masses using QCD sum 
rules. 
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